solve_triangular_system Interface

public interface solve_triangular_system

Module Procedures

private pure function solve_triangular_system_mtx(a, b, upper) result(rst)

Solves a triangular system of the form where is a triangular matrix, either upper or lower, for equation .

Arguments

Type IntentOptional Attributes Name
real(kind=real64), intent(in), dimension(:,:) :: a

The N-by-N triangular matrix.
real(kind=real64), intent(in), dimension(:,:) :: b

The N-by-NRHS matrix.
logical, intent(in), optional :: upper

An optional argument specifying if the matrix is upper triangular (true), or lower triangular (false). The default assumption is that is an upper triangular matrix.

Return Value real(kind=real64), allocatable, dimension(:,:)

The N-by-NRHS solution matrix, .

private pure function solve_triangular_system_vec(a, b, upper) result(rst)

Solves a triangular system of the form where is a triangular matrix, either upper or lower, for equation .

Arguments

Type IntentOptional Attributes Name
real(kind=real64), intent(in), dimension(:,:) :: a

The N-by-N triangular matrix.
real(kind=real64), intent(in), dimension(:) :: b

The N-element array.
logical, intent(in), optional :: upper

An optional argument specifying if the matrix is upper triangular (true), or lower triangular (false). The default assumption is that is an upper triangular matrix.

Return Value real(kind=real64), allocatable, dimension(:)

The N-element solution array, .